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ABSTRACT 

We report on optical imaging of the X-ray binary SAX J1808.4— 3658 with the 8-m Gemini South 
Telescope. The binary, containing an accretion-powered millisecond pulsar, appears to have a large 
periodic modulation in its quiescent optical emission. In order to clarify the origin of this modulation, 
we obtained three time-resolved r'-band light curves (LCs) of the source in five days. The LCs can 
be described by a sinusoid, and the long time-span between them allows us to determine optical 
period P — 7251.9 s and phase 0.671 at MJD 54599.0 (TDB; phase 0.0 corresponds to the ascending 
node of the pulsar orbit), with uncertainties of 2.8 s and 0.008 (90% confidence), respectively. This 
periodicity is highly consistent with the X-ray orbital ephemeris. By considering this consistency and 
the sinusoidal shape of the LCs, we rule out the possibility of the modulation arising from the accretion 
disk. Our study supports the previous suggestion that the X-ray pulsar becomes rotationally powered 
in quiescence, with its energy output irradiating the companion star, causing the optical modulation. 
While it has also been suggested that the accretion disk would be evaporated by the pulsar, we argue 
that the disk exists and gives rise to the persistent optical emission. The existence of the disk can be 
verified by long-term, multi-wavelength optical monitoring of the source in quiescence, as an increasing 
flux and spectral changes from the source would be expected based on the standard disk instability 
model. 

Subject headings: binaries: close — stars: individual (SAX J1808.4— 3658) — X-rays: binary — stars: 
low-mass — stars: neutron 



1. INTRODUCTION 

While it had long been believed that neutron 
star (NS) X-ray binaries (XRBs) are progen- 
itors of the recycled millis e cond radio pulsars 
(|Bhattacharya fe van den Heuvell 119911 ). it was the 
discovery of coherent pulsations from the transient 
XRB SAX J1808.4-3658 (hereafter J1808.4) dur- 
ing its X-ray outburst in 1998 that first and finally 
confirmed the connection betw een the two systems 
(jWiinands fe van der Kl is 1998): in this binary, the 
accreting NS is a 2.49 ms X-ray pulsar. As the first 
example of accretion-powered millisecond pulsar sys- 
tems, J1808.4 has been extensively studied, with various 
interesting properties revealed (see lHartman et aT1l2008l 
and references therein). In this paper, we focus on the 
optical periodic modulation seen in this binary and 
report on our observational study of the modulation. 

The orbital period of J1808.4 is P orb -7249.157 
s (~2.01 hr), accurately known to one part in 10 10 
from Do ppler modulations of the millisecond pul- 
satio ns (fChakraba rtv fc Morgan! 119981 : lHartman et al.l 
2008(1 . Combined with the derived mass function 
of 3.8xlO _5 M0, the period implies that the mass- 
transferring companion could be a 0.17 Mq low-mass 
main-s equence star, but more likely a ~0.05 Mq brown 
dwarf (Bildstcn & Cha krabartvll2001 ) . At a distance of 
D = 3.5 kpc (|Gallowav fc Cumm ing 2006|), the optical 
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counterpart in quiescence is several magnitudes brighter 
(V = 20.7, Ly ~ 3.0 x 10 32 ergs s _1 assuming isotropic 
emission and extrinction Ay = 0.73; see § 2 and § 4) 
than the possible types of stars suggested as the compan- 
ion, probably indicating that the op tical emission arise s 
from the accretion disk in the binary ([Homer et al.l2001ft . 
However in the quiescent state, 10-40% sinusoidal-like 
modulations in the source's optical light curves (LCs) 
have been reported (jHomer et al.ll200lt ICampana et alJ 
2004), and this is puzzling because the quiescent X- 
ray lumi nosity is approxima tely Lx — 5 x 10 31 ergs 
s _1 (e.g., iHeinke et al.l l2007h . two orders of magnitude 
lower than that requ ired to account for the modulation 
(|Burderi et al.ll2003t) . Typically in a low-mass X-ray bi- 
nary (LMXB), sinusoidal optical modulation arises from 
X-ray heating of the companion star by the central X- 
ray source: the visible area of t he heated face varies as a 
function of orbital phase (e.g., lArons fc King||T993[ ). In 
J1808.4, depending on the companion's star types, only 
0.5-1.4% [estimated by (R 2 / D b ) 2 / A, where i? 2 is radius 
of the companion and is the separation distance be- 
tween the NS and companion] of the total energy flux 
from the central NS would be received by the compan- 
ion for isotropic emission. This leads to the suggestion 
that in quiescence, the NS might switch to be a rotation- 
powered pulsar so that the rotational energy would be 
the energy source that heats half su rface of the cqmpan - 
ion star and causes the modulation (|Burderi et al.f 2003) . 

However, there are other possibilities that do not re- 
quire a rotation-powered pulsar, and we have consid- 
ered whether or not the accretion disk could give rise to 
the modulation. It has been known that "superhumps" , 
which are com monly seen in short-period cataclysmic 
variables (CV; IWarnerl 1 1995T ) , also appear in LMXBs 
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(e.g.. lHaswell et aT1l2001l) . These periodic modulations 
have periods a few percent longer than the orbital periods 
and can be sinusoidal-like with an amplitude of ~10%, 
arising from a processing , eccentric accretion disk (e.g., 
IWhitehurst k K ing 1991f l. Indeed, it has been suggested 
that those NS LMX Bs with P nrh < 4.2 h r are potential 
superhump sources ([Haswell et all 120011) . In addition, 
several parts of an accretion disk co uld contribute sig- 
nifica ntly to optical modulation fe.g. JMason k Cordoval 
1982). It has also been suggested that for an X-ray tran- 
sient, its quiescent optical emissi on may come from a 
bright spot on the accretion disk (jMenou fc McClintockl 
120011) . 

In particular, the superhump possibility was suggested 
by the X-ray LC obtained in the source's 2002 outburst. 
As shown in Figure [TJ the LC exhibits a ~5-day peri- 
odic modulation at the end of the outburst. If this in- 
dicates the precession periodicity (-P prC c — 5 days) of 
the accretion disk, it would imply a superhump period 
of P sh = 7373 s (1/P sh = 1/Porb - 1/Pproc) and super- 
hump excess e = 0.017 [e = {P s h — Pmb) / Poib] (jPattersonl 
2001). The excess value is consistent with those obtained 
for c a taclysmic variables and LMXBs ([Patterson et al.l 
2005; Ha swell et al.l 1200 ll) . Furthermore, a mass ra- 
tio of q ~ 0.08 cou ld be estimated from the relation 
e = 0.18(? + 0.29g 2 ([Patterson et al.ll2005h . implying a 
companion mass of 0.11 M for 1.4 Mq NS mass. This 
companion mass is within the range implied by the mass 
function. 

Previously, time-resolved imaging observations over a 
small period of time (covering only ^1.5 orbital pe- 
riod of the binary) were made. However, these obser- 
vations were carrie d out either with a small telescope 
([Homer et al.l I2001D or under very poor observing con- 
ditions ( Campana et al.l|2004f ). resulting in large uncer- 
tainties in the obtained LCs. In order to study the op- 
tical emission from J1808.4, and particularly to probe 
whether it could be a superhump source, we have ob- 
tained high quality optical LCs of the source in its qui- 
escent state through time-resolved photometry. The ob- 
servations were made with the 8-m Gemini South Tele- 
scope over five days, allowing us to determine the period 
and phase of the o ptical modulati on accurately. We note 
that lHeinke et alJ (|2008f ) (see also lDelove et al.ll2008l ) re- 
cently observed the source simultaneously at X-rays and 
optical g'i' wavelengths, and from the observations they 
confirmed the inconsistency between the large amplitude 
optical modulation and low X-ray luminosity. 

2. OBSERVATIONS AND DATA REDUCTION 

To determine the periodicity in the source's optical 
emission accurately, three Gemini queue mode observa- 
tions of J1808.4 were carried out in five days, on 2008 
May 11, 12, and 15. The starting time of each observa- 
tion was approximately 06 hour (UTC) each day, result- 
ing in a time span of ~4 days between the first and third 
observations. We proposed such observations because we 
estimated that the time span would allow us to determine 
the period to <10 s accuracy, and the second observation 
would be needed to keep the track of the optical period- 
icity phase. A Sloan r' filter, with the central wavelength 
at 6300 A, was used for imaging. The detec tor was the 
Gem ini Multi-Object Spectrograph (GMOS; iHook et all 
l200l . which consists of three 2048x4608 EEV CCDs. 
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Fig. 1. — X-ray light curve of J1808.4 during its 2002 outburst, 
obtained with the Proportional Counter Array (PCA; 2-60 keV 
energy range) on board the Rossi X-ray Timing Explorer satellite. 
Gaps in the light curve were due to Earth occultations of the source. 
At the end of the outburst, a ~5 day periodicity is tentatively 
suggested. 

We used 2x2 binning, providing a pixel scale is 0.146" 
pixel -1 . In each night, 36 images of the source were 
obtained contiguously, each with approximately 3.9 min 
exposure time. The detector's slow read mode, having 
55 s readout time, was used. As a result, the total obser- 
vation time in each night was approximately 3 hrs, cov- 
ering 1.5 orbital cycles of J1808.4. The average seeing 
[full-width half-maximum (FWHM) of the point spread 
function (PSF) of the images] for the three nights were 
0.63", 0.58", and 0.70", respectively. The second night 
had the best observing conditions, with the seeing reach- 
ing 0.51" a few times during the observation. 

We used the Gemini IRAF package GMOS for data 
reduction. The images were bias subtracted and flat 
fielded. The bias and flat frames were from GMOS base- 
line calibrations, made on 2008 May 13 and May 11, 
respectively. The st andard star used f or flux calibration 
was PG1047+003A ([Smith et al.l[2002] ). The observation 
of this star was made on 2008 May 13, also as part of the 
GMOS baseline calibrations. The airmass of the obser- 
vation was 1.234, which can be estimated to have caused 
a zero-magnitude offset of 0.03 mag 4 . We did not add 
this offset to our brightness measurements given below; 
instead we consider it as an uncertainty for absolute flux 
calibration. 

We performed PSF-fitting photometry to measure the 
brightnesses of the source a nd other in-field stars . A 
photometry program D0PH0T ([Schechter et al.lll993h was 
used. A finding chart of the target is shown in Figure O 
As can be seen, our target is located between two stars 
with similar brightnesses. Its distance to star a is 0.6" 
and to star b is 1.0". In a few of images, we have FWHM 
around 0.8"; in these cases, our target and star a are 
nearly unseparated. For these images, we positionally 
calibrated them to a reference image that was combined 
from four best-quality images in night 2. We determined 
the positions of our target and star a in the reference 
image and fixed them at the positions for photometry of 
the images. 

We performed differential photometry to eliminate sys- 
tematic flux variations in the images. An ensemble of 8 
isolated, nonvariable stars in the field were used. The 

4 www.gemini.edu/sciops/instruments/gmos/calibration/photometric- 
stds 
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Fig. 2.— Gemini South r' image of the J1808.4 field. Object X, 
located between star a and b, is the optical counterpart to J1808.4. 
The star labeled as C is used as a check star. 

brightnesses of our target and other stars in each image 
were calculated relative to the total counts of these stars. 
Star C (Figure^]) was used as a check star, because it was 
non- variable and had similar brightness to our target. 

We used the third image from the second night to ob- 
tain absolute magnitudes of the target and nearby stars, 
as it is one of the best-quality images. The aperture cor- 
rection was calculated using 15 in-field stars, with an un- 
certainty of 0.025 mag. The resulting magnitudes of the 
target are given in Table [TJ and the average magnitudes 
of the nearby stars a and b, and the check star C were 
r' = 21.492 ± 0.048, 21.133±0.013, and 21.178±0.013 
mag, respectively. The LCs of our target and stars a and 
C are shown in Figure[3l As can be seen, star a was likely 
a variable, with the magnitudes and standard deviations 
of its three LCs being 21.545 ± 0.029, 21.444 ± 0.020, 
and 21.486 ± 0.021. The difference between the first and 
second nights is 2.9cr significant. These results are sum- 
marized in Table [2] 

As we compared our results with those previously re- 
ported, we noted that the source magnitudes, resulting 
from imaging observations made on 1999 July 11 with 
the Very Large Telescope (VLT) at the European South- 
ern Observatory, are approximate l y 1 m ag lower than 
the values given bv lCampana et all (|2004), who analyzed 
the same data. The data consist of 1 min exposures in 
the V, R, and /-bands, taken with the high resolution 
collimator, providing a pixel scale of 0.1" pixel -1 . The 
instrumental magnitudes were calibrated against photo - 
metric standard stars in the SA110 field ()LandoliJll992f ). 
We obtained V = 20.73 ± 0.04, R = 20.59 ± 0.04, and 
I = 20.15 ± 0.06, where the uncertainty is the quadratic 
sum of the uncertainty in the zeropoint, the aperture 
correction, and the instrumental magnitude. Comparing 
the magnitudes of the in-field stars, including star a and 
6, from the VLT observations and ours, we believe that 
our magnitude values are correct. 

3. PERIODICITY DETERMINATION 

As can be seen in Figure [31 the LCs of J1808.4 
clearly show a sinusoidal modulation, and appear to 



have different average brightnesses, indicating overall 
variations from day to day. The times of the data 
points are barycentric corrected, with the JPL Solar 
System Ephemeris DE405 used. In order to determine 
the modulation, we fit the LCs with function m — 
m c +mh sin[27r(t/P+</>o)]j where t is the time, P, <j>o, and 
rah are the period, starting phase, and semiamplitude of 
the sinusoidal modulation, respectively. The parameters 
m c and rrih were kept as a constant for each LC, but were 
allowed to have different values in different LCs. As a re- 
sult, we found that the best-fit sinusoid (x 2 — 1879 for 
100 degrees of freedom) has P = 7251.9 s and 4>q = 0.671 
at MJD 54599.0 (TDB; Phase = 0.0 corresponds to the 
ascending node of the pulsar orbit). 

While the LCs can be described by the sinusoidal func- 
tion, as shown in Figure [3j the large x 2 value indicates 
large scattering of the data points from the best-fit func- 
tion. There is a systematic uncertainty caused by our 
target's proximity to star a. This can be seen from the 
fact that the standard deviations of the three LCs of 
star a are significantly larger than its uncertainties from 
PSF-fitting (the average is 0.013 mag) and the standard 
deviation (0.013 mag) of all data points of the check star 
C. In addition, we also independently used the program 
DA0PH0T in the ESO-MIDAS system for photometry. The 
resulting LCs are very similar to those resulting from 
D0PH0T, but with the standard deviations of the differ- 
ences between the two sets of the LCs being 0.027, 0.019, 
and 0.014 mag for the three nights. These values are ap- 
proximately consistent with the standard deviation val- 
ues of star a, confirming the contamination of the pho- 
tometry caused by the proximity of our target and star 
a. Adding the standard deviations of star a in quadra- 
ture with the uncertainties of data points (resulting from 
PSF-fitting) of the target, the x 2 value is reduced to 266 
for 100 degrees of freedom. This indicates that there is 
intrinsic scattering of the data points from the single si- 
nusoid. For example, we note that the brightest data 
point in each LC appears at phase 0.05-0.17 after the 
maximum of the sinusoid. This pattern is likely to be 
true, because the D0PH0T and DA0PH0T measurements at 
the LCs' region are nearly identical. 

The uncertainty on P is 2.8 s (90% confidence), found 
from Monte Carlo simulations. We generated 10,000 sets 
of simulated LCs, each like the sets of the actual data 
points. In doing that, we used the best-fit parameters 
and added to each set of LCs Gaussian-distributed de- 
viates, where the Gaussian distribution was estimated 
from the residuals to the best-fit model. Having standard 
deviation a — 0.04 mag, the Gaussian mimics the rela- 
tively large scattering of the data points from the best-fit 
model. We then fit each set of simulated LCs with a si- 
nusoidal function. The uncertainty on P was determined 
by the spread of values. We also determined the uncer- 
tainty on the phase this way, and found it to be 0.008 
(90% confidence). Comparing to the X-ray ephemeris 
(phase at MJD 54599.0 is ~0 6 714 with a negligible un- 
certainty; |Hartman~£F!aIl2008|) , the optical periodicity is 
consistent with being orbital. We investigated whether 
the period uncertainty might be caused by the uncer- 
tainty on the GMOS exposure recording, because it is 
not clear how accurate the latter was. We made simu- 
lations by assigning randomly produced, uniformly dis- 
tributed time offsets to the recorded image times, and 



4 



found that the period value is not sensitive to any pos- 
sible offsets. For example, conservatively assuming 1-s 
uniformly distributed offsets for the GMOS time record- 
ing, the resulting period difference has a range of 0.03 s, 
negligible compared to the statistical period uncertainty. 

The average brightness of J1808.4 in the three nights 
increased from 2f.l23, to 21.105, to 21.023 mag, while 
the semiamplitude of the modulation decreased from 
0.214, to 0.202, to 0.191 mag (Table©. These variations 
may suggest that the two components of the emission, 
the persistent and modulated, were independent of each 
other; as the former was increasing, the modulation frac- 
tion was decreasing. However the uncertainties on these 
parameters are relatively large, ~0.04 mag (90% confi- 
dence), showing that the variations of the semiamplitude 
are not significant. This is because each of our obser- 
vations covered only 1.5 orbital cycles, insufficient for 
an accurate determination. Therefore we conclude that 
we have detected an approximately 20% flux modulation 
from J1808.4 in r' band. In addition, the optical peaks 
correspond to when the pulsar is right in front of the 
companion (superior conjunction of the companion; 270° 
mea n orbital longi t ude), confi rming the previo us results 
from lHomer et all (|2001h and lCampana et all (|2004h . 

4. DISCUSSION 

Using the 8-m Gemini South Telescope, we have ob- 
tained, for the first time, well-determined LCs from 
J1808.4 in its quiescent state over a time span of four 
days. From the above studies of the LCs, we find that 
the optical period and phase are consistent with the X- 
ray ephemeris, indicating that the optical modulation is 
orbital in origin. In studie s of several tens of LMXBs 
at op tical wavelengths (e.g.. Ivan Paradiis fe McClintockl 
1995), in no instance has there been an accretion disk 
giving rise to a sinusoidal modulation at the orbital pe- 
riod. In addition, the sinusoidal maximum must cor- 
respond to superior conjunction of the companion star. 
Because of these, we rule out the possible disk origin for 
the modulation that we have suspected. However, the 
source in outburst could still be a superhumper, which 
might have been hinted in the X-ray LC (Figure [l}. As 
the outward extension of accretion disks in outburst has 
both been o bserved and reproduced in disk instability 
simulations jOsaki 1993; iDubus et al.l l200lh - it would 
not be unexpected for the accretion disk in J1808.4 to 
have extended to the resonance zone during the 2002 
outburst, develo ping into an eccentric for m due to the 
tidal instability (|Whitehurst fe Kind[l99lD . In fact, su- 
perhumps have been seen i n outbursts of both black-hole 
and NS LMXB sys tems (jO'Donoghue fc Charles! Il99fil : 
lElebert et al.l [2008) . In order to determine this pos- 
sibility for the long periodicity seen in the 2002 out- 
burst, time-resolved imaging observations, like ours, of 
the source in outburst are needed. Since the source will 
be as bright as ~17 mag in an outburst fe.g-. lGiles et al.1 
Il999t IWang et al.ll20 0lf). a search for superhump modu- 
lation will be feasible even with a small telescope. 

Bas ed on the current observational stu dies of LMXBs 
(e.g., Ivan Paradijs fc McClintockl Il995f ) . it seems ex- 
tremely unlikely that the observed optical modulation 
would arise from a source other than the compan- 
ion star. Thus far, pulsar wind heating of the com- 



panion is the only model that has bee n suggested 
(|Burderi et all 120031 : ICampana et~aT1 l2004f ). The long 
term spin-down rate of the pulsar has been measured, 
indic ating a rotational en ergy loss rate of 9xl0 33 ergs 
s" 1 (jHartman et al.l |2008T >. This energy output, pre- 
sumably in the form of a pulsar wind, would illumi- 
nate the companion star. Assu ming isotropic emission 
and a brown dwarf companion (|Bildsten fc Chakrabartvl 
|2001[ ). the fraction of the total energy received by the 
companion is ~0.005?7»(i?2/0.13 i?©) 2 , where 77, is the 
fraction of the r eceived energy absorb ed by the compan- 
ion. Following lArons fc King! (|1993h . the companion's 

heated face would have temperature ~7430??y 4 K, due 
to pulsar wind heating by the putative rotation-powered 
pulsar. Using such a hot face that varies following a 
function of [1 + sini sin(27rf/P)], where i is the incli- 
nation angle of the binary, and also including a con- 
stant flux component Fc, we tested whether we could 
re-generate the averaged LCs of J1808.4. The distance 
and extinction to the source were fixed at 3.5 kpc and 
Ay = 0.73, respectively, where the extinction value is 
estimated from Ay — 7Vh/0.179 x 10 22 cm -2 by as- 
suming hydrogen column density to the source N-g = 
0.13 x 10 22 cm" 2 (jDickev fcLockmanll 19901: iHeinke etldl 
20071). The extinct ion law for Sloan filters given by 
Schlegel et all ([19981) was used. We found that the pa- 
rameter values of i ~ 63° (M 2 ~ O.O49M ), 77* ~ 0.46, 
and Fc — 19 /xJy can provide the observed modulation 
(the resulting \ 2 — 2100, with no systematic uncertain- 
ties considered). Although we used a very simple model, 
these derived parameter values are consistent with its 
known properties. In addition to the fact that the com- 
panion is likely a ~0.05 M Q star, the source shows no 
X-ray eclipses or dips, implying i < 70°. The obtained 
77* values ar e within the range found for two binary ra - 
dio pulsars (|Stappers et al.ll200lURevnolds et al.ll2007f) . 
in which it is known that the companion is irradiated by 
the pulsar wind. Therefore, it is plausible that the NS in 
J1808.4 does turn into a rotation-powered pulsar in qui- 
escence, giving ris e to the optical mo dulation. We note 
that very recently, iDeloye et al.l (|2008l ) used an advanced 
model to fit their g'i' light curves, and also found that 
the required heating energy should be ~ 10 34 ergs s _1 , 
consistent with the d erived spin-down lum inosity (which 
has 30% uncertainty: lHartman et al"1l2008l ). 

The origin of th e pers istent optical emission is not 
clear. iHomer et al.l (|2001[ ) tried explaining the emission 
from an X-ray irradiated disk around the pulsar, but it 
may not be appropriate to use a steady thin disk model 
to describe a disk in the thermally stable cold state (lower 
cold branch of the standard thermal equilibrium S'-curve; 
e.g., lLasotal 12001] ). since a disk temperature profile in 
the cold state can be drasti cally different from th e hot 
state (the steady disk case). ICampana et al.l (|2004l ) used 
a shock front, arising from the interaction between the 
companion star and pulsar wind, and the irradiated com- 
panion to account for the emission. Here we argue that 
the accretion disk in quiescence exists, against the sug- 
gestion that the disk would be evapora ted by the pulsar 
(Burden et aLll2003t iHeinke et al.ll2008h . and this can be 
tested by monitoring J1808.4 at optical wavelengths. 

Accordin g to th e standard disk instability model (DIM; 
e.g., lOsakil Il996t lLasotal I200l1) . while the mass accre- 
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Fig. 3. — Optical r' light curves of J1808.4 (diamonds), in which sinusoidal modulation is clearly visible. For comparison, the LCs of 
the nearby star a (triangles) and check star C (circles, downward shifted by 0.8 mag) are also shown. The best-fit sinusoidal function to 
the LCs of J1808.4 is shown as the solid curves, while the constant magnitude for each LC is indicated by the das hed lines. The optica l 
periodicity well matches the X-ray ephemeris (dotted curves), which gives the mean orbital longitude of the binary (Hartman ct al. 2008). 
The optical brightness peaks correspond to when the pulsar is right in front of the companion star (270° mean orbital longitude). The 
brightest data points in the LCs, indicated by arrows, are at phase 0.05-0.17 after the maxima of the sinusoid. 



tion rate to the NS in J1808.4 is very low during qui- 
escence, M acc < 6.2 x 10 -15 Mq yr _1 (estimated from 
the observed X-ray flux), the average mass transfer rate 
from the companion to the accretion disk is as high as 
^10 _11 Mq yr -1 (estimated fr om the X-ray fluence in 
each outburst; iGallowavl 120081 ) . The transferred mass 
is stored in the disk, building up the surface density 
for triggering the next outburst. The average persis- 
tent r' flux from J1808.4 in our observations is esti- 
mated to be Fc = 19 /iJy corresponding to a disk 
luminosity of L r > — 2ttD 2 Fc/ cosi ~ 1.5 x 10 32 ergs 
s _1 (i = 63° is assumed). There is plenty of gravita- 
tional energy available to power this emission as matter 
moves inwards through the outer disk. At the time av- 
erage accretion rate of M ~ 10~ n M Q yr _1 , matter 
falling into a radius of 4000 km releases gravitational 
energy at a rate that matches the observed luminos- 
ity. This radius is far larger than those that are sug- 
gested for the inner radius r ln of the disk. Generally, 
ri n would be close to the Alfven radius, r; n ~ 56 km 
(M in /lO- u M yr- 1 )- 2 / 7 (Ai/lO 26 G cm 3 ) 4 / 7 , where M in 
is the mass accretion rate in the inner edge of the disk and 
H is the magnetic moment, /i ~ 10 26 G cm 3 for J1808.4 
(jHartman et al.1 l2008h . In the cold state, M- ln would 
be lower than M, and we note that for M ln — 0.1M 
(|Dubus et al.ll200l[ ). r; n is 110 km. However, since a ra- 
dio pulsar presumably would have no interactions with 
a surrounding disk, r- m would be larger than the light 
cylinder radius of the pulsar, which is 120 km. As can be 
seen, it is possible that in quiescence, the disk in J1808.4 



would be outside of the light cyliner. In addition, the 
disk temperature profile in quiescence may be des cribed 
by a c onstant, at least rig ht after an outburst (e.g.. lOsa"kl 
119961 iDubus et al] 120011 ). For J1808.4, we find that an 
effective temperature of 4600 K for the disk can give rise 
to the persistent r' flux, where the disk is assumed to be 
cut off at the tidal radius 3.7xl0 10 cm (~ 0.9i?i, where 
i?i is the NS's Roche lobe radius). This temperature 
value is consistent with those typically considered in the 
DIM (Lasota 2001; the critical effective temperature for 
having an outburst is ~ 6000 K). 

In order to verify our suggestion that the persistent 
optical emission arises from the disk, long-term, multi- 
wavelength optical monitoring of the source in its qui- 
escent period is required. From such observations, we 
might expect to see an increasing flux from the source. 
Moreover, since in the DIM the temperature profile 
as a function of disk radius is predicted to be chang- 
ing, turning from a constant right after an outburst 
to a power-law-like func tion prior to the next outburst 
(e.g., IDubus et al.1l2001| ). we would also see flux spec- 
trum changes. This type of well-behaved changes would 
not be expected from the pulsar wind shock model 
(|Campana et al.l l2004f ) . thus allowing to determine the 
origin of the persistent emission. 

If the companion star is irradiated by the pulsar wind, 
there is no reason to think that the disk is not. It 
has been suggested that the disk in quiescence might 
be evaporated by the pulsar (e.g., iBurderi et al.l [20031 : 
iHeinke et al.l I2008D. but according to the recent calcula- 
tions by Jonesl (|2007l ). a pulsar wind may only be effective 



() 



in heating a disk. Basically, as X-rays from a NS would 
ionize the surface of a disk, the Poynting flux, which is 
dominant in a wind when it is not far from the light cylin- 
der of the pulsar, would interact with the ionized parti- 
cles, conver ting energy i nto disk heating. Using equa- 
tion (16) in Uonesl (|2007f ). we estimate that the baryon 
loss rate of the disk at the inner radius is approximately 
3xl0 21 (r in /120 km)" 3 citT 2 s" 1 , only 0.05% of the sur- 
face density (~10-100 g cm -2 ) that is g enerally con- 
sider ed in the accretion disk models (e.g., iDubus et al.l 
120011) . This suggests that the disk in J1808.4 could exist 
and might be irradiated by the pu l sar w ind. However, 
using the model provided by Uonesl (|2008l ). the flux due 
to pulsar wind heating would be 2 /xJy for parameter 
C = 0.3 (0.03< C ^0.3 and a larger £ value correspo nds to 
a hig her disk effective temperature; see details in Uonesl 
120081) . The flux would be 10% of the average r' flux, 
which would suggest a weak pulsar-wind heating effect 
in J1808.4. 

Finally, it will be of great interest if J1808.4 can be 
determined to become rotation-powered during quies- 
cence. We note th at the source could be very similar to 
PSR J2051-0827 (IStappers et al.lll996f ). a binary mil- 
lisecond pulsar system. For example, the latter has an 
orbital period of 2.38 hr and a mass function of l.OxlO -5 
Mq , and the pulsar has a spin-down luminosity of 6 x 10 33 
ergs s _1 . However, searches for pulsed r adio emission 
from J1808.4 have not been successful (e.g. jBurgav et all 



2003). Here we suggest that the source might be identi- 
fied by searching for its pulsed 7-ray emission. Observa- 
tions of millisecond pulsars suggest that t heir efficiency 
at 7- ray energies may be as high as ~7% (|Kuiper et al.l 
2000). This implies a 7-ray flux of ~5xl0 -13 ergs cm -2 
s _1 for J1808.4, possibly detectable by deep observations 
with Fermi Gamma-Ray Space Telescope. 
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TABLE 1 
Photometry of J1808.4 



MJD a r' Ar' b 



0.257907 21.028 0.006 

0.261323 21.116 0.006 

0.264715 21.069 0.006 

0.268084 21.148 0.009 

0.271469 21.267 0.010 



Note. — Table 1 is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. 
a Days since MJD 54597. 0. b la uncertainty resulting from 
PSF fitting. 



TABLE 2 

Summary of brightnesses of nearby stars and J1808.4 

IN OUR OBSERVATIONS 



Source 


Obs 1 


Obs 2 


Obs 3 




(MJD 54597) 


(MJD 54598) 


(MJD 54601) 


Star a 


21.545 ±0.029 


21.444 ±0.020 


21.486 ±0.021 


Star b a 




21.133±0.013 




Star C* a 




21.178±0.013 




Sinusoidal fitting 


J1808.4 








Average magnitude 13 


21.12 


21.11 


21.02 


Scmiamplitude b 


0.21 


0.20 


0.19 



Note. — Uncertainties of 0.025 mag and probable 0.03 
mag from the aperture correction and zero point calibration, 
respectively, are not included. 

a Average magnitude is derived from all three 
observations. 13 Uncertainties (90% confidence) are ~0.04 
mag. 



